Specific heat evidence for two-gap superconductivity in ternary-iron silicide Lu2Fe3Si5 
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We report low-temperature specific lieat studies on single-crystalline ternary-iron silicide super- 
conductor Lu2Fe3Si5 with Tc = 6.1 K down to '~rc/20. We confirm a reduced normalized jump in 
specific heat at Tc, and find that the specific heat divided by temperature C/T shows sudden drop at 
•^Tc/5 and goes to zero with further decreasing temperature. These results indicate the presence of 
two distinct superconducting gaps in Lu2Fe3Si5, similar to a typical two-gap superconductor MgB2. 
We also report Hall coefficients, band structure calculation, and the anisotropy of upper critical 
fields for Lu2Fe3Si5, which support the anisotropic multiband nature and reinforce the existence of 
two superconducting gaps in Lu2Fe3Si5. 

PACS numbers: 74.25.Bt, 74.70.Dd, 71.18.-|-y 



Since the discovery of high- Tc superconducting inter- 
metallic compound MgB2 with Tc ~ 39 K [l|, the mate- 
rials containing Ught elements, such as, B, C, and Si, has 
attracted much attention. Among compounds containing 
light elements, the ternary-iron silicide superconductors 
i?2Fe3Si5 with R = Lu, Y, Sc, Tm, or Er, have stimulated 
many investigations of the relation between superconduc- 
tivity and magnetism due to magnetic atoms 0, S 
i?2Fe3Si5 crystallizes in the tetragonal Sc2Fe3Si5-type 
structure consisting of quasi-one-dimensional iron chain 
along c^axis and quasi-two-dimensional iron square struc- 
tures [5;]. The properties of these compounds are very 
peculiar in several aspects. For instance, one of the high- 
est superconducting transition temperature among any 
known iron compounds are observed in Lu2Fe3Si5 with 
Tc = 6.1 K and Sc2Fe3Si5 with T, = 4.5 K 0,0. In addi- 
tion, a large negative pressure effect on Tc is reported in 
Lu2Fe3Si5 and Sc2Fe3Si5 (dTc/dp = -7 x 10"^ K/ bar) 
while a large enhancement of Tc is exhibited in Y2Fe3Si5 
{Tc/dp = -33 X 10"^ K/ bar) f^. Moreover, reentrant 
superconductivity is found for Tm2Fe3Si5 [3] and the co- 
existence between superconductivity and antiferromag- 
netism is observed in Er2Fe3Si5 

One of the most puzzling issues in the ternary-iron 
silicide superconductors is the peculiar superconducting 
state in Lu2Fe3Si5. In this material, significant non- 
BCS features, such as a large residual linear term in 
the superconducting specific heat and a reduced normal- 
ized specific heat jump at Tc smaller than BCS value of 
1.43, have been reported [1, SUBl- In addition, a non- 
exponential power-law specific heat below Tc has been 
observed y,y,|lo|- Moreover, a rapid decrease in Tc when 
a small amount of nonmagnetic impurities replace some 
of Fe atoms [llj . These anomalous features let us specu- 
late a possibility of the spin-triplet superconductivity in 
Lu2Fe3Si5. However, the Josephson effect measurement 
between Lu2Fe3Si5 and Nb is against the triplet pairing 
state [13]. It has been an unsettled problem whether 
unusual superconducting properties of Lu2Fe3Si5 are in- 



trinsic or arise from magnetic impurities due to struc- 
tural disorder and defects. In this letter, we address 
the issue of the unusual superconducting properties of 
Lu2Fe3Si5 by the low temperature specific heat measure- 
ments of high-quality single crystals. We report that the 
peculiar temperature dependence of specific heat can be 
well-described by the phcnomcnological two-gap model 
similar to a typical two-gap superconductor MgB2. We 
also investigate the Hall effect, band structure, and the 
anisotropy of upper critical fields for Lu2Fe3Si5, in or- 
der to support the validity of two-gap model in terms of 
electronic structure. 

Single crystals are grown by the floating-zone tech- 
nique using an image furnace in Ar atmosphere at a rate 
of ~2 mm/h. The starting polycrystalline rod is prepared 
by fusing several polycrystalline lumps, prepared by arc 
melting the constituent elements in Ar atmosphere. The 
single crystals are annealed for 4 weeks at 800° C fol- 
lowed by 5 days at 1250°C. It should be noted that the 
annealing at high temperature for an extended period of 
time is essential to obtain sharp superconducting tran- 
sition. Magnetization was measured by a commercial 
SQUID magnetometer (MPMS-XL5, Quantum Design). 
The Hall coefficient measurements are performed by stan- 
dard six-wire configuration. Specific heat was measured 
by the relaxation method in a ^He refrigerator. A band 
structure calculation is carried out by using a full po- 
tential APW (FLAPW) method with the local density 
approximation (LDA) for the exchange correlation po- 
tential. 

Figure [1] depicts the specific heat divided by temper- 
ature C/T as a function of T^ for Lu2Fe3Si5. A jump 
in the specific heat at 6.1 K indicates that the transition 
into the superconducting state is very sharp, which is also 
observed in magnetization measurements shown in the in- 
set of Fig. [TJ Above Tc, C/T slightly deviates from T^- 
dependence. The normal-state specific heat, which keeps 
the entropy balance at Tc, is fitted by using the equation, 
C = jnT -\- PnT^ + anT^ , whcrc 7„T is the electronic 
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FIG. 1: (color online) The specific heat divided by temper- 
ature C/T as a function of for Lu2Fe3Si5. The dashed 
hne represent the fit to the data in the normal state C — 
JnT + pnT^ + a„T^. The soUd line is the fit to the data 
in the superconducting state C = 7sT + PsT^ ■ Inset shows 
temperature dependence of the zero-field-cooled (ZFC) and 
field-cooled (FC) magnetization for Lu2Fe3Si5 at 5 Oe. 
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FIG. 2: (color online) Temperature dependence of the nor- 
malized electronic specific heat for Lu2Fe3Si5. Dashed and 
solid lines are BCS normalized specific heat and a two-gap 
fit calculated by Eqs. ([1} and ([2]), respectively. Inset shows 
semi-logarithmic plot of C^./'y-nTc vs Tc/T for Lu2Fe3Si5. The 
solid line is a linear fit below ~ Tc/5. 



term and PnT^ + q;„T'^ represents the phonon contribu- 
tion. We obtain the parameters 7„ — 23.7 mJ/mol K^, 



origm. 

We now discuss the detailed temperature dependence 
of specific heat for Lu2Fe3Si5. Figure [2] shows temper- 
ature dependence of the normalized electronic specific 
heat for Lu2Fe3Si5. For comparison, we also plot the 
BCS weak-coupling limit specific heat (dashed line) cal- 



culated based on the following equations 14| . 
S 6 Ao 



/3„ = 0.276 mJ/mol K^, and a„ = 1. 73x10^3 mJ/mol ^„T, 
K^, which are very close to the values for the previous 
polycrystalline samples [1, 0, [l^ . The large T^-term ob- 
served in the normal-state specific heat suggests com- 
plex phonon density of states and has been reported for 
other ternary superconductors The magnitude of 

the jump at Tc, AC, is 150 mJ/mol K. The value of the 
normalized specific heat jump at I^C j'^nTc-, is 1.05, 
significantly smaller than the BCS value of 1.43, close to 
reported values d, [^, [l3] ■ Below the specific heat for 
Lu2Fe3Si5 shows peculiar T-dependence. At high tem- 
peratures from ~ Tc to ~ Tc/b, C/T shows almost per- 
fectly r^-dependence. We estimated a linear extrapola- 
tion term, 7s, from Tc to T = by using the equation, 
C = -fsT + l3sT^. From the fit to the data plotted in Fig. 
[11 we obtain the parameters 7^ — 12.5 mJ/mol K'^ and 
(is = 1.37 mJ/mol K^. The linear term 7s is about 53 
% of the normal-state electronic coefficient. These val- 
ues are again consistent with those reported previously 
B S [13] ■ The observation of the reduced normalized 
jump at Tc and the large linear extrapolation term in 
the superconducting specific heat for both polycrystalline 
and high-quality single crystalline samples strongly sug- 
gest that these anomalous features stems from intrinsic 
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where / = (exp(/3£;) + l)'^ and /3 = {keT)-'^. The en- 
ergy of quasi-particles is given hy E — {e^ + A^(t))^/^, 
where e is the energy of the normal electrons relative 
to the Fermi energy, and the temperature dependence 
of an energy gap varies as A(t) — AoS{t), where S{t) 
is the normalized BCS gap at the reduced temperature 
t = T/Tc as tabulated in Ref.[l3|. The integration vari- 
able is y = e/Ao. Apparently, temperature dependence 
of Ce/jnT for Lu2Fe3Si5 is very different from that of 
BCS specific heat. In the BCS weak-coupling limit, below 
Tc, Ce/7„T drastically decreases with decreasing tem- 
perature and is almost zero below Tc/5. In contrast to 
the BCS weak-coupling fimit, Ce/7„T for Lu2Fe3Si5 de- 
creases slowly with decreasing temperature from Tc to 
~ Tc/5. Ce remains about 60 % of normal-state specific 
heat even at ~ Tc/5. Below ~Tc/5, CehnT shows a sud- 
den second drop and goes to zero with further decreasing 
temperature. Very similar behavior is also observed in 



MgB2 [IJ, [16|, llT], ll8|. In order to explain the anoma- 



lous behavior of specific heat for MgB2, two-gap model 
has been proposed and it reproduces the peculiar spe- 
cific heat well 3, l3]- We show here similar analysis for 
Lu2Fe3Si5. Within the two-gap model, the total specific 
heat can be given by the sum of the contribution from 
each gap calculated independently using the Eqs. ([T]) and 
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FIG. 3: (color online) Temperature dependence of the Hall 
coefficients for Lu2Fe3Si5 with (a)/ || a,H \\ c and (b)/ || 
c, H II a. 



(12). The solid line shown in Fig. [2] represents the two-gap 
fit with the magnitude of two distinct superconducting 
gaps, 2Ai/kBTc = 4.4 and 2/^2/kBTc = 1.1, and the ra- 
tio of weight, 71 : 72 = 47 : 53, where a partial electronic 
specific heat coefficient 7^ with 7n = 71 -I- 72 character- 
izes each band. The calculation based on the two-gap 
model shown in Fig. [2] well reproduces the experimental 
results, supporting the presence of distinct two super- 
conducting gaps in Lu2Fe3Si5 similar to MgB2. At high 
temperatures just below Tc, the magnitude of AC/7„r is 
determined by 71 because quasi-particles excited across 
the smaller gap A2 have negligible contribution to the 
specific heat jump at T^. By contrast, at low tempera- 
tures T ^ Tc, C/T shows the sudden decrease and its 
magnitude is determined by 72 because the smaller gap 
A2 determines the excitation of quasi-particles. There- 
fore, all the anomalous features observed in the supercon- 
ducting specific heat for Lu2Fe3Si5, such as the reduced 
jump of the specific heat at Tc, the large residual linear 
term, and the sudden decrease below '■^ Tc/b, can be well 
described by the two-gap model. Similar results are ob- 
tained from other single crystalline samples. It should be 
noted that penetration depth measurements also support 
the two-gap superconductivity in Lu2Fe3Si5 19]. 

A different approach reinforces the existence of smaller 
superconducting gap. In the weak-coupling BCS super- 
conductors, it is well known that empirical relation Ce oc 
exp(-aTc/T) with a 1.44 holds for 2.5 < Tc/T < 6 
[20| . The inset of Fig. [2] shows a semi-logarithmic plot 
of Ce/"fnTc vs Tc/T. Solid line represents a linear fit 
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FIG. 4: Fermi surfaces of Lu2Fe3Si5 calculated by a FLAPW 
method: hole-like Fermi surfaces from (a) 155th and (b) 156th 
bands, and (c) electron-like Fermi surface from 157th band. 



to the data at low temperatures following the relation 
Ce/lnTc ~ exp(-aTc/T) with a = 0.29. Obviously, the 
observed value of a is about four times smaller, suggest- 
ing that a fully opened gap exists at low temperatures 
but its magnitude is very small. All these features are 
observed in 

The temperature dependence of Hall coefficients also 
supports that Lu2Fe3Si5 is a multiband system. Fig- 
ures [3l[a) and (b) depict the temperature dependence 
of Hall coefficients Rh for Lu2Fe3Si5. In the configu- 
ration with I \\ a and || c, at high temperatures near 
300 K, the sign of Rh is negative. Rh begins to in- 
crease with decreasing temperature from ~150 K. At ^ 
110 K, the sign of Rh changes to positive. With fur- 
ther decreasing temperature, Rh decreases after showing 
a maximum at around 40 K. At low temperatures, Rh 
increases again with decreasing temperature after show- 
ing a minimum. In the configuration with / || c and 
H \\ a, the sign of Rh is negative in the whole tempera- 
ture range. At high temperatures, Rh decreases with de- 
creasing temperature and the slop becomes slightly larger 
below ~ 150 K. Rh becomes nearly constant below 70 
K and increases steeply with decreasing temperature be- 
low 20 K. In a simple metal, the Hall coefficient in sin- 
gle band system is written as, Rh = ±l/ne, where n 
is a carrier number. Rh is almost T- independent. In 
contrast, the Hall coefficient in multiband system, for in- 
stance, consisting of electron and hole bands, is given by, 
Rh = ("/iMh - / {e{nhtXh + ?^eMe)^), where Uh (ne) 
is a carrier number of hole (electron) band and /i/i (/Ze) 
is a mobility of hole (electron) band. The Hall coefficient 
in multiband system can be temperature dependent and 
even show the sign change because of the different tem- 
perature dependence of the mobility in each band. The 
strong temperature dependence of the Hall coefficient for 
Lu2Fe3Si5 in both configurations strongly suggests that 
Lu2Fe3Si5 is a multiband system. 

The Fermi surfaces in MgB2, where two-gap supercon- 
ductivity has been studied extensively, is formed by two 
kinds of bands: two-dimensional cr-band of B and three- 
dimensional TT-band of B [2l|, [23, [2^ • It is believed that 
the larger gap opens in cr-band while the smaller gap 
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opens in 7r-band. Figure |4] shows the Fermi surfaces for 
Lu2Fe3Si5 calculated by an FLAPW method. The Fermi 
surfaces consist of two hole-like bands (155th and 156th) 
and one electron-like band (157th), each having contribu- 
tion to the density of states at the Fermi level iV(0), 13.9 
%, 44.5 %, and 41.6 %, respectively. In the two hole- 
like bands, some part of Fermi surfaces are quasi-one- 
dimensional, reflecting the iron-chain structure, while 
others are three dimensional. Unfortunately, it is not 
clear which parts of the bands shown in Fig|4] are re- 
sponsible for the larger and the smaller gaps. From the 
analogy of MgB2, however, we speculate that the exis- 
tence of anisotropic band and three-dimensional one may 
play an important role in the two-gap superconductiv- 
ity. Evaluation of Fermi surface dependent gap as has 
been done for MgB2 is necessary to clarify this point 
24| . From A^(0) obtained by the band calculation, we 
can obtain the calculated specific heat coefficient jtand 
= 8.69 mJ/mol K^, which is about three times smaller 
than 7„ in the present study, indicating that presence of 
the strong electron-phonon and/or electron-electron in- 
teractions. 

Figure [5] shows temperature dependence of upper crit- 
ical field along a-axis and c-axis obtained by resistivity 
measurements. The temperature dependences of Hc2 in 
both directions are almost linear in the present tempera- 
ture region. Apparently, as shown in Fig. [5l H^2 is larger 
than H^2 i^i the whole temperature range. Anisotropy 
of upper critical field H^2/Hc2 is ~ 2, indicating that 
Lu2Fe3Si5 is a weakly one-dimensional superconductor 
reflecting the crystal structure. In fact, the zero-field 
inter-plane resistivity (/ || c) is about four times smaller 
than in-plane resistivity (/ || a) just above Tc as shown 
in the inset of Fig. O 

In summary, we have prepared high-quality single crys- 
tals of ternary iron-silicide superconductor Lu2Fe3Si5 
and performed the low-temperature specific heat studies 
down to ~ Tc/20. The measurements reveal a reduced 
normalized specific heat jump at Tc, a large linear term 
even at ^ Tc/5, and a second drop of Ce/T below Tc/5. 
These results can be well reproduced by the phenomeno- 
logical two-gap model similar to a typical two-gap su- 
perconductor MgB2. The Hall coefficient measurements 
and band structure calculations support the multi-gap 
superconductivity. We emphasize that Lu2Fe3Si5 gives 
a unique opportunity to investigate the multi-gap super- 
conducting properties in the following reasons: (1) sim- 
plicity of the multi-gap nature, i.e., each band has nearly 
the same contribution to density of states, (2) relatively 
high transition temperature, and (3) large crystal size. 
These advantage, as compared with several other possi- 
ble multigap superconductors, such as MgB2, NbSe2 [1^, 
PrOs4Sbi2 l26|], URu2Si2 ^ may shed hght on the de- 
tails of multi-gap superconductivity. 
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FIG. 5: (color online) Temperature dependence of upper crit- 
ical field fioHc2 along a-axis (•) and c-axis (o) obtained by 
resistivity measurements. Dashed lines are linear fits to the 
data. Inset shows the temperature dependence of zero-field 
inter- (7 || c) and in- plane (/ || a) resistivity. 
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